Boundary effects on one-particle spectra of Luttinger liquids 
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We calculate one-particle spectra for a variety of models of Luttinger liquids with open boundary 
conditions. For the repulsive Hubbard model the spectral weight close to the boundary is enhanced 
in a large energy range around the chemical potential. A power law suppression, previously predicted 
by bosonization, only occurs after a crossover at energies very close to the chemical potential. Our 
comparison with exact spectra shows that the effects of boundaries can partly be understood within 
the Hartree-Fock approximation. 



Interacting fermions in one spatial dimension do not 
obey Fermi liquid theory. It took about thirty years 
to fully understand the generic low-energy physics of 
one-dimensional (ID) fermions with repulsive interac- 
tion, now called Luttinger liquid (LL) behavior ||l[. For 
various correlation functions LL theory predicts anoma- 
lous power laws with interaction dependent exponents [Q . 
Several systems were studied using different experimen- 
tal techniques in order to verify these predictions. Pho- 
toemission or transport measurements seem to be most 
promising ||l|-^. While the basic understanding of LL 
behavior emerged in the study of ID systems with peri- 
odic boundary conditions (PBC), it is obvious that in or- 
der to describe realistic geometries more realistic bound- 
ary conditions have to be used. Therefore hard walls, 
usually called "open" boundary conditions (OBC) were 
studied recently For the translational invariant 

system standard renormalization group (RG) arguments 
can be used to show that for repulsive interactions 
the 2fci?-scattering part (usually called "gi-interaction" ) 
scales to zero and the Tomonaga-Luttinger model (TLM) 
describes the generic low energy physics. In some of the 
previous publications on the open boundary problem [^|J^] 
it was tacitly assumed that similar RG arguments can also 
be applied to the system without translational invariance. 
Therefore a quadratic form in boson operators was used 
to describe the electron-electron interaction. Then it is 
straightforward to calculate correlation functions 1^,^ . It 
was shown that the local spectral density p{x,u}) enter- 
ing the description of angular integrated photoemission 
is strongly modified near the endpoints of a ID chain. 
The algebraic vanishing of the spectral density with fre- 
quency to close to the chemical potential jj, was found 
to be governed by a coupling dependent boundary expo- 
nent ub which, for repulsive interaction, is larger than 
the bulk exponent |^,^. This result was used recently in 
the interpretation of photoemission experiments |^,^ . In 
Ref. as has been studied for a variety of integrable 
models using results of the Bethe ansatz (BA) and as- 
suming conformal invariance. 



Thus far only specialized methods (bosonization, the 
BA and conformal field theory) have been used to study 
the one-particle properties of LL's with boundaries and 
assumptions were necessary to justify the applicability of 
these methods. Furthermore these methods do only pro- 
vide the exponent which characterizes p(x, ui) at asymp- 
totically small energies. In this letter we study several 
models for LL's with OBC using bosonization and the 
density matrix renormalization group (DMRG) method. 
By comparison with results obtained from these (nu- 
merically) exact methods we will show that the low en- 
ergy behavior of p{x, to) can partly be understood within 
the (non-self-consistent) Hartree-Fock (HF) approxima- 
tion for the self-energy, which already gives power law 
behavior. This is a very surprising observation as the 
HF approximation for the case of PBC does not capture 
any of the LL features. For the models considered we 
will explicitly demonstrate that as can be expressed in 
terms of the bulk LL parameter Kp, as — i^p^ ~ 
where z = 1 in the spinless case and z = 2 for spin 1/2- 
fermions Q-^. For the TLM and the lattice model of 
spinless fermions (SF) with nearest neighbor interaction 
the exponent found in HF agrees to leading order 
in the coupling with as and for small couplings there is 
quantitative agreement between the exact local Green's 
function and the HF approximation. We will show that 
for the Hubbard model (HM) with small or moderate 
Coulomb repulsion U the spectral weight at the endpoint 
of the chain is enhanced in a large energy range close to 
fi. A crossover to a power law decrease occurs in a sur- 
prisingly small energy range around fi. This crossover 
behavior cannot be obtained from bosonization or confor- 
mal field theory due to the limitation to asymptotically 
small energies. Because of the "flowing gi-coupling" 0| 
there is only qualitative agreement between DMRG and 
the HF approximation in case of the HM. Besides their 
relevance for LL's with OBC our results are important 
for the investigation of LL's with PBC including impu- 
rities, since such models are expected to scale to chains 
with OBC g. 
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We first present a careful discussion of the interaction 
term of the Hamiltonian in the basis of the one-particle 
eigenstates (pnix) = ^2/Lsin (fc„a;) with fc„ — nn/L, 
n S IN for OBC, a continuum model, and an interaction 
with a spatial range R. Here L denotes the length of 
the chain. Note that the kn do not have the meaning of 
momenta. In terms of the creation and annihilation oper- 
ators of the eigenstates the two-body interaction 
takes the usual form |^ 

n,n' ,m,m' 

with the matrix elements Vnmn'm' ■ If we express products 
of the sine-functions tpn in terms of cosine-functions they 
read 

F {kn k-a' 7 k-fa -|- fcjn' ) F (/c^ -\- kn' , k^i ^m') 
-\-F {kn ~\~ kn' , kni ~\~ kni' 

)]/L, (1) 

with 

F{q,q') = Y f dx f dx' cos{qx)V{x — x')cos{q'x'). 
L Jo Jo 

Due to the fact that kn = mr/L (and not 2'Kn/L) in 
addition to terms involving Kronecker deltas a finite size 
correction appears for non-local interactions 

F{q. q') = V{q)/2 {6,^^, + Sg,-g,) + g{q, q')/L, (2) 

where V{q) is the Fourier transform of the interaction 
V{x). The correction term g{q, q')/L does not contribute 
to a_B and we thus neglect this term. For an in- 
teraction which is long range in real space, considered 
by Tomonaga for PBC, V{k) is different from zero for 
I A; I < fcc, where the cut-off kc = 1/i? is much smaller 
than the Fermi momentum kp = npi^/L. Then only 
the first term on the right hand side (rhs) of Eq. (|l|) 
is important for the low energy physics and V can be 
written as a bilinear form in boson operators. Using 
bosonization correlation functions can then be calculated 
as discussed in Ref. Q yielding a_B = {Kp^ — l)/z with 

Kp = [l + zV{Q)/{-kvf)]~''^'^ and the Fermi velocity vp- 
The full uj dependence of p{x,uj) for the TLM has been 
presented in Fig. 1 of Ref. Q. 

In contrast to the bulk exponent, which, for weak cou- 
pling, is quadratic in the interaction, as has a contri- 
bution linear in the interaction. Thus signs of the non- 
analytic behavior of p{x, w) can already be obtained using 
the HF self-energy. The delta terms on the rhs of Eq. (||) 
yield 

(^^^)^^,^Skk'lsp^^ J2 [v{k~k,) + V{k + k^) 

I ki<kF 



+ 1. |^y(fc + k')^V {[k - k']/2)} f {[k + k']/2) (3) 

{^^(^ -k')-V {[k + k']/2)] f {\k - fc'l/2) , 

where /(fcm) is the Fermi function, which is equal to one 
if m is integer and m < np and zero otherwise, and Sfi = 
zV{Q)nF/ L. The calculation of p^^ requires a matrix 
inversion. The leading non-analytic behavior in the HF 
approximation can be traced to the term in which is 
proportional to /([fc -I- k']/2). To leading order in V the 
term in the second line of Eq. gives a contribution 
to proportional to [V{0) - zV{2kF)] log |a; - 
For a long range repulsive interaction V{2kp) = and 
the prefactor of the logarithm is positive. This indicates 
the suppression of spectral weight close to /i^^. Note 
that the logarithmic divergence is not due to a singu- 
lar frequency behavior of E^^, as it occurs in second 
order perturbation theory for PBC, but emerges in the 
process of the matrix inversion. For finite systems it is 
numerically straightforward to perform the matrix inver- 
sion and calculate p^^. Instead of comparing p^^ and 
the exact spectral function p obtained from bosoniza- 
tion as a function of to we have studied the spectral 
weight at the chemical potential and position x, denoted 
wo{x,nF',V), for a given kp as a function of n^?. The 
ratio wo{x, np; V)/wo{x, np; 0) as a function of 1/np dis- 
plays the same kind of power law behavior as does p{x, uj) 
as a function of uj. For the lattice models considered 
later we will only be able to exactly determine the spec- 
tral weight at the chemical potential on a lattice site 
I as a function of the system size. In Fig. ^ we com- 
pare the 1/np dependence of wo{x,np;V)/wo{x,np;0) 
for the exact solution and the HF approximation. Al- 
ready the HF approximation displays power law behav- 
ior, which shows that the leading logarithmic divergence 
of p^^ discussed above can be resummed to give a power 
law. A detailed study shows that = y(0)/(27njHF) 
where uhf = vp + y(0)/(27r) is the Fermi velocity in 
the HF approximation and thus and as do agree 
up to leading order in ^(0)7(27^^). Quantitative agree- 
ment between exact results and HF can be reached for 

V{0)/{27TVp) < 1. 

For the lattice models (SF, HM) considered next even 
the case with OBC can be solved exactly by the BA, but 
as for PBC not much about correlation functions can be 
learned directly from the BA. Information about bound- 
ary exponents can be obtained if conformal invariance is 
assumed Q . First we will discuss the local spectral func- 
tion at site i = 1 in a lattice model of SF with lattice 
constant a = 1, hopping matrix element t = 1, nearest 
neighbor interaction U, and OBC. We have performed a 
DMRG study @ for chains of up to iV = 512 sites calcu- 
lating matrix elements wo{np;U) = \{Eq''~^\ci\E^'')\'^ , 
i. e. the spectral weight at the chemical potential and the 
boundary site. Here I-Eq^) denotes the exact ni?-particle 
groundstate and ci destroys a fermion at site 1. Typi- 
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cal data for wo{nF',U)/wo(nF]0) at quarter filling and 
for U = 1 are shown in Fig. |l|. For these parameters 
Kp — 0.8447 as can be found in Ref. |11|. For large np 
the numerical data nicely follow the solid line, which is 
proportional to a power law with exponent as — 0.1838 
and we can thus conclude, that for SF the suppression is 
described hy as- 
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FIG. 1. Spectral weight at /i and close to the boundary. 
The circles show the exact results for the spinless TLM at 
X = 3R with a long range interaction V'(g) — UQ{kc — \q\), 
U/{2ttvf) = 0.1 and kp = 4fcc. The dashed line presents the 
HF approximation. The squares show the exact results for 
the lattice model of SF at site i = 1 for np/N — 0.25 and 
[7 = 1. The respective HF approximation is given by the 
long dashed line. The solid line displays a power law with 
exponent ap = 0.1838. The error of the DMRG data is of 
the order of the symbol size. 



The HF self-energy for SF is best studied in the site 
representation. Due to the OBC the site occupancies 
{rii) and the renormalization of the hopping amplitudes 
{cl_^_lCi) depend on i and in order to obtain the site 
diagonal Green's function one numerically has to in- 
vert a tridiagonal matrix. Fig. |^ presents HF data for 
wo(nF]U)/wo{nF; 0) and the same parameters as above. 
The HF approximation shows power law behavior with 
exponent = T4fi/(27rwHF) which has the same form 
as for the TLM if one replaces V{0) by the effective inter- 
action Foff = V{0) - V{2kF) = 2U[1 - cos(2fci.)]. Again 
aB and do agree up to order Vcs/{2ttvf) and both 
curves show quantitative agreement for Vcs/i^nvp) 1. 

Next we will consider the HM and first discuss the re- 
sults of the Hartree approximation (the Fock term van- 
ishes). To obtain the site diagonal Green's function one 
again has to invert a tridiagonal matrix which we have 
done numerically for systems of up to = 10^ lattice 
sites. The corresponding spectral function pi{uj) at the 
boundary site i = 1 is shown in Fig. ^ for various values 
of the onsite interaction U and fixed filling np/N = 0.4. 
We have connected the individual weights of the finite 
system to a continuous line. In contrast to the expecta- 



tion from the general statement in previous work the 
spectral weight near ji^ is strongly enhanced. This could 
have been expected already from Eq. (^). The prefactor 
F(0) — zV{2kp) of the leading logarithmic correction to 
^HF -g gjvgj^ i-,y (I _ z)V{Q) for a fc-independent delta 
interaction, i.e. in a model with spin it has the opposite 
sign as in the case of a long range interaction. As long 
as \U / {2-Kv p)]\og {N) ^ 1 this indicates a logarithmic 
increase of the weight close to ji^ . If we go to larger 
values of [C//(27rwi?)] log (A^), i.e. larger A^ (or [/), as in 
Fig. ^, we observe a crossover to a power law decay for 
energies extremely close to ji^ . It can be shown that 
the crossover occurs at energies which are exponentially 
(in —\/U) close to /i^. In Fig. || the crossover cannot 
be observed as the energy resolution Aw ~ 1/A^ is too 
low. In the inset of Fig. || the spectral weight of the oc- 
cupied states is shown on a log-log scale. The increase 
of weight close to /i^ is given by a power law with expo- 
nent —U I {2ttvp) which is cut off at the crossover energy. 
This shows that the logarithmic divergence obtained in 
leading order in U can be resummed to produce a power 
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FIG. 2. Spectral density at the boundary site for the HM 
with N — 800 sites and np — 320 in the Hartree approxima- 
tion and for different values of U. 




FIG. 3. The same as in Fig. 2, but for [/ = 3, iV = 10^ 
and np = 4 - 10^. The inset shows the spectral weight for the 
occupied states uj < on a log- log scale (solid line). The 
dashed lines display power laws with exponents ±U / {2-nvp). 
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law in a large energy range close to /z^. The expo- 
nent of the subsequent power law suppression is given 
by Ul{2nvF). 

In order to demonstrate that the crossover behavior 
found in the Hartree approximation is present also in 
the exact spectral function we again have calculated ma- 
trix elements 'Wo{np;U) using DMRG for chains of up 
to = 256 sites. In Fig. ^WQ{np;U)/wo{np;0) is pre- 
sented for np/N = 0.25, U = 0.5,2.5 and [/ = 8. For 
small N and U the DMRG and the Hartree results do ap- 
proach each other as expected. A detailed study shows 
that the difference between the weights for fixed N goes 
like U"^. If one goes from small to large N the Hartree 
result for U = 2.5 diplays the power law increase with 
exponent —U/{2ttvp), the crossover, and the subsequent 
decrease with exponent U/{2'KVp) similar to pi{to) as a 
function of w. For U — 0.5 only the increase can be 
seen as the crossover occurs at much longer chain length 
not shown in Fig. IJ. The exact weight obtained from the 
DMRG also increases for increasing N. From the U = 2.5 
data it is clear that the crossover occurs at smaller N and 
the maximum is less pronounced compared to the Hartree 
result. For U = 8 the weight already decreases even for 
the smallest system sizes available. The dashed-dotted 
line displays a power law with exponent as- The bulk 
LL parameter Kp for a given U and np/N can be found 
in Ref. The DMRG data are consistent with the 

conclusion that the final power law suppression near the 
boundary is given by the boundary exponent as 
The increase of weight in the U = 0.5 DMRG data fol- 
lows a line which is proportional to a power law with 
exponent —as- The leading behavior of the boundary 
exponent is as ~ U/{Attvp) p3| which is one-half oi the 
exponent U/{2'KVp) found in the Hartree approximation. 
This kind of discrepancy between exponents obtained in 
perturbation theory and the leading behavior of exact 
exponents is known from the HM with PBC and occurs 
because a fixed point is reached only when backscattering 
has scaled to zero 0. For U ~ 2.5 the crossover length 
is too small to observe a well defined power law increase. 
We can conlcude that analogous to the Hartree solution 
the DMRG results are consistent with an increase of spec- 
tral weight for small [U/{2Trvp)] \og{N), a crossover and 
a subsequent power law decrease. 

In conclusion for the models studied we have explic- 
itly confirmed the expression for the boundary exponent 
as in terms of the bulk LL parameter Kp obtained from 
bosonization and the BA. This exponent characterizes 
the one-particle spectral weight for asymptotically small 
energies and close to the boundary. Investigating the 
spectral function at higher energies we have found a quite 
unexpected behavior for the HM. For small and moderate 
positive U the spectral weight close to the boundary first 
increases when the energy approaches the chemical po- 
tential from below and only in a narrow energy range the 
power law decrease with exponent as sets in. Whether 



LL behavior has convincingly been demonstrated in pho- 
toemission experiments is still a matter of debate 
Nonetheless the above scenario has to be taken into con- 
sideration in any discussion of the influence which bound- 
aries (open ends or impurities) might have on the pho- 
toemission spectra of quasi-one-dimensional conductors. 
We have furthermore demonstrated that the influence of 
boundaries on one-particle spectra of LL's can partly be 
understood within the HF approximation. 
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FIG. 4. Spectral weight at /i and the boundary site for 
the quarter filled HM. The symbols show the results obtained 
from DMRG. The dashed curves display the weight in the 
Hartree approximation. The solid line shows a power law 
proportional to {1/N)~°'^ for U = 0.5 and the dashed-dotted 
line is proportional to {1/N)°'^ for U — 8. The error of the 
DMRG data is of the order of the symbol size. 
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